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We study Andreev bound states in a long shallow normal constriction, which is open to 
a superconductor at both ends. The interesting features of such setup include the absence 
of electron-hole symmetry and the interference of electron waves along the constriction. We 
compare results of a numerical approach based on the Bogoliubov equations with those 
of a refined semiclassical description. Three types of Andreev bound states occur in the 
constriction: i) one where both electron and hole wave part of the bound state propagate 
through the constriction, ii) one where neither electron nor hole wave part propagate, and 
Hi) one where only the electron wave propagates. We show that in a wide energy region the 
spacing between the Andreev states is strongly modulated by the interference of electron 
waves passing the constriction. 



I. INTRODUCTION 

Andreev bound states^ occur in a normal metal system of finite size that is in contact with bulk superconductors. 
They exhibit complicated and sometimes counter-intuitive physics which depends much on disorder and electron 
scattering in the normal metal. The most adequate method of describing Andreev states in metals is a semiclassical 
one.El In many occasions the method can be reduced to a simple circuit theory.0 

A new challenge is presented by the fabrication of systems where the normal part is represented by a doped 
semiconductor. Typically, Andreev states are being formed in a two-dimensional electron gas. There are means to 
shape such a 2DEG in different ways forming controllabe constrictions. There are also means to provide conditions 
for ballistic transport of the electrons, so that they are scattered at the constriction boundaries-only. The Andreev 
states in such systems can also be addressed by the traditional semiclassical technique, see e.g.Q However there are 
cases where this approach fails. 

The point is that the traditional semiclassical approach essentially exploits an approximate electron-hole symmetry. 
This symmetry holds provided the superconducting energy gap A exceeds a typical energy of the electrons. Re- 
cently it has been pointed out that this condition may fail in few-mode superconductor-semiconductor-superconductor 
constrictions. LJ 

If electron-hole symmetry does not hold, the problem should be adressed by means of the Bogoliubov equations 
that provide a strict quantummcchanical treatment. The equations can be easily solved for a separable geometry, 
That is why the occurrence and character of Andreev bound states in a "straight" SNS junction are well understood.^ 
A two-dimensional picture of such a system is shown in FIG. |l|a. The simplicity of the problem comes from the simple 
shape of the boundaries. 

As soon as the boundaries are curved, the geometry is not separable any more. The quantummechanical problem 
appears to become hardly solvable acquiring features of quantum chaos. In the present paper we present a simple 
model system where complications related to chaos may be lifted but the absence of the electron-hole symmetry plays 
an important role. We have chosen for a SNS junction with a constriction in the middle, schematically depicted 
m FIG. 0b. We assume the adiabaticity of the confining potential that allows to reduce the problem to a set of 
one-dimensional problems for non-mixing transverse modesE We solve the one-dimensional Bogoliubov equations 
numerically. To understand the results, and the degree of their generality, we present in addition two semiclassical 
methods for the one-dimensional problem. 

The paper is organized as follows. In Section || we derive the Bogoliubov equations for the adiabatic constriction. 



We specify the model in use and analyze related energy scales in Sec. III. We present a simple semiclassical picture 
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of Andreev levels in the Sec. IV. A brief 



ualitative description of the results is gi ven i n Sec. The numerical 
We give concluding remarks in Sec. VI] . An interpretation of the 
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results are presented and discussed in Sec. 

interference features found is given in Appendix [X], along with the underlying semiclassical description 
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FIG. 1. A straight SNS junction and a SNS constriction. 



II. THE BOGOLIUBOV EQUATIONS FOR THE MODEL CONSTRICTION 



The Bogoliubov equations for a clean system with free electrons in the normal metallic part of the system can be 



written as 
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The two-component wave function ^(r) describes quasiparticle excitations of the superconducting state, the energy E 
of those is being counted from the Fermi energy \x. The spatial dependence of the gap function A(r) remindes that we 
deal with an inhomogeneous system. In the superconductors A(r) is assumed to be a complex constant, while in the 
normal part of the system A(r) = 0. We restrict ourselves to effectively two-dimensional systems. This is achieved 

technically by choosing the thickness w in one direction, let us say the z direction, such that ^(^) 2 > A 4 , by which 
transverse modes in the z direction are forbidden. _ 

In the following derivation we go along the path proposed irJH for normal constrictions. We assume that the boundary 
potential in the remaining transverse direction, the y-direction, is infinitely high, so that the wave function must be 
equal to zero at the transverse boundaries, which can be expressed by a sine function. The total wave function 'J, 
which is now a function of two coordinates, then can be presented as the linear combination of the following form 
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(j) n being a two-component wavefunction. Substituting this form in Eq. (|l|) and taking the inner product with 
sin ^^fy(y — ^2^)) > on e obtains the equations for the upper component of 4> n 
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The corresponding equation for the lower component differs from Eq. Q3J) only by a minus sign of the left hand 
side. The components are decoupled, since in the constriction the superconducting gap parameter A equals zero. For 
the time being we will work with a slowly varying transverse shape of the constriction, and assume that the terms 
containing derivatives d (x) and d (x) can be neglected. After we have specified the system we will give an estimate of 
these terms. In this approximation, the different transverse modes are uncoupled. Therefore we obtain the following 
equation for each mode 
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FIG. 2. The shape of the model constriction that corresponds to a linear potential U(x), for two different Lt 



where we have defined an effective potential U(x) = ^(^jj) > representing the transverse mode energy. The two- 
component wave function in the longitudinal direction of the system, 4> n (x), assumes different forms in the different 
parts of the system. In the superconductor left of the normal metal it has the form 
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|Al| and k 2 Fmax = 2m^/h 2 — ^ 77 , d maa; being the width in the superconductor. Similarly, at the 
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right hand side of the constriction this form becomes 
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Obviously, the labels L and R refer to the location of the superconducting part with respect to the constriction. 
The above forms shall be used as boundary conditions for Eq. (El). 



III. MODEL 



The model constriction is depicted in FIG. |^. The constriction is symmetric. To simplify numerical calculations, 
we choose the shape of the boundaries in such a way that the effective one-dimensional potential U (x) either remains 
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FIG. 3. The reflection coefficient from the top as function of energy. 



constant of varies linearily. This choice has the advantage that the solutions of Eq. (Q) can be readily expressed in 
terms of well-known Airy functions. The effective potential and the corresponding shape of the constriction, d(x), are 
shown in FIG. |[ 

The figure shows only half of the junction, the part in the negative x direction. The constriction is long so that 
Lkp ^S> 1. The constriction is shallow so that the top of the potential U max for a transverse mode of interest almost 
matches Fermi energy, Umax ~ fJ-- We thus define SE = fi — U m ax, \SE\ <C fi. This is achieved by a proper choice of 
the width d t in the middle of the constriction. An important parameter is also L t , the length of the narrowest part 
of the constriction. The width reaches its maximum value d max at |x| = (L — L t )/2, which is at a distance L s from 
the superconductors. 

Despite apparent simplicity of the model, it exhibits a variety of relevant energy scales. We list them going from 
bigger ones to smaller ones. The biggest scale is obviously the Fermi energy itself, \i. Another important scale we 
obtain by considering electron reflection from one side of the top and disregarding existence of another side. Classically, 
the electron is fully transmitted at energies exceeding U max and fully reflected otherwise. Due to quantum effects, 
there is a smooth change from full reflection to full transmission that occurs within an energy interval E c . There is 
a related length scale x c that corresponds to minimal electron wavelength. We estimate the typical length scale x c 
and the energy scale E c by equating potential and kinetic energy: E c ~ (dU/dx)x c ~ h 2 /mx 2 . Let us note that by 
virtue of our model the potential derivative can be estimated as dU/dx ~ n/L. By that we define an energy scale 
E c = 0.03(j,(kpL)~ 2 / 3 at which the reflection of electron waves from the top is changing to full transmission. The 
numerical coefficient is choosen in such a way that the reflection probability reaches 0.5 at E — U max « E c . The 
energy dependence of the reflection probability is plotted in FIG. The curve roughly behaves as R ~ E c /E. 

Three different energy scales arise from the interference of the electron waves either beyond or at the top of the 
potential. Beyond the top, the typical wavevector is of the order of kp and the interference condition reads 6kL ~ n. 
This gives En ~ Ti 2 kpSk/m TTfi(kpL) . The typical traversal time for this part of the system is of the order 



of h/En. The wavevector at the top is considerably smaller, fc 



top (mSE) 2 /h. This may lead to much longer 
E c we obtain the corresponding energy scale 



traversal times through the top of the constriction. Substituting SE 
E l2 = {2L/L t )(E c /^) 1 / 2 E ll , and E i2 « E n if L ~ L t . 

The third energy scale we obtain by estimating the energy counted from the top at which a few wavelengths match 
the length of the top, kt op Lt — f. This yields Eis ~ n(kpL t )~ 2 . The same energy scale E^ determines when 
tunneling through the top becomes essential. .-.We see that Ei2 -C E c , at least if L ~ L t . Usual assumptions for 
adiabatic constrictions are such that Ei2 ~ E C B This makes our theory distinct from, for instance, RefB. 

The above analysis of our model captures the essential features of electron propagation in long shallow junctions 
and all qualitative results are actually model- independent. 

We have performed numerical calculations for constrictions that have a length L ranging from 10 5 k] 
and d max = at the SN boundaries. The width at the middle d t is chosen to vary between 17 anc 

assures that U max ~ fi for the transverse mode n = 6. 

Using the Airy function solutions of Eq. (Q), we estimate the errors involved in neglecting the terms in Eq. (||), 
which contain derivatives of d(x). As far as the first derivative is concerned, looking at the second and third term in 
Eq. (H), we calculated the relative errors 
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for the two Airy functions Ai(x) and Bi{x), for three x values, at the beginning of the constriction, somewhere in the 
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FIG. 4. Three types of semiclassical Andreev states. The type III is specific for the absence of electron-hole symmetry. 



middle, and at the position where the narrowing has come to an end. Looking at the fourth term, and realizing that 
for our linear potential U(x) — tx + c the relation d d = 3d holds, we calculated the relative error 
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at the same positions as well. For the three points mentioned and for five different energies smaller than |A| the 
average error in the first derivative appeared to be smaller than 0.01%, and the error in the second derivative was 
even a factor of 10 3 smaller. It could be concluded, that the omitted terms are really negligible. 



IV. SEMICLASSICAL ANDREEV STATES 



Let us start by presenting a semiclassical picture of Andreev states in the simplest case where electrons and holes are 
either fully transmitted or fully reflected from the potential U(x). We count energies of the states from the chemical 
potential. We denote SE = fi — U max . 

We assume no phase difference between superconductors. The Andreev state with positive energy E is made of 
electron waves with energies E and — E, those we call "electrons" and "holes" respectively. There may be three 
situations, depicted in FIGjJ: I. Both electron and hole are fully transmitted through the constriction, SE > and 
E < 5E. II. Both electron and hole are reflected from the top, SE < and E < \5E\. III. The electron is transmitted 
and the hole is reflected, E > \SE\. Since the potential U(x) varies slowly, the quantized energies of Andreev states 
may be obtained from the Bohr-Sommerfeld quantization rule applied to Eq. (||). The solutions of Eq. (0) can be 
written in the form 

e * and e ifk.( s )d s (g) 

for the electron and hole part of <fi n (x) respectively, where k± = yj2m{p — U(x) ± E)/h. We begin our discussion with 
the case I. In this case, an electron travelling from left to the right builds up a phase of f_ L / 2 k+(x)dx. The electron 
undergoes Andreev reflection and returns as a hole, and the wave function gains the additional phase cf>A — acos(E/A), 
see Eqs. (||) and (||). The hole builds up a phase of J L ^ 2 k-(x)dx. After Andreev reflection at the left side, the 
particle has made a complete roundtrip. The quantization condition is that the phase gain is a mupliple of 2tt, that is 

,L/2 

/ (fc+(x) - k-{x))dx + 24> A = 2nn (10) 

J-L/2 
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Differentiating Eq. (10) with respect to E gives an approximate relation for the energy difference between adjacent 
Andreev states, 



E n+1 - E n = 2irh{T el + T hole + 2h-^y 1 (11) 

oE 

where we have introduced semiclassical times of travelling from the one interface to the other one, for electrons and 
holes respectively, 

f L ' 2 dx f L ' 2 dx , , 

T e i = / ; and T ho i c = / ; (12) 

' J-L/2 ^2(M - U(x) + E)/m J-L/2 y/2(» -U(x)-E)/m V ' 

and made use of semiclassical relation between the travelling time and the derivative of the quantummechanical 
phase with respect to the energy. If there are many bound Andreev states, their typical spacing shall be small being 
compared to A, so that T e i t hole 3> K/A. Since the derivative of the Andreev phase is of the order of 1/A, the second 
term in the right hand side of Eq. (|ll| ) can be disregarded and E n+ i — E n as 2irh(T e i +T/ lo ; e ) _1 This suggests a 
presentation of our numerical results: we will plot the energy spacing times the travelling time as a function of the 
state energy. For the sake of symmetry we will use 

^{E n+1 -E n ) Av{T ol + T holc } = 6, (13) 

where Av{T cl + r hole } = \{T c \(E n ) + T hole (E n ) + T cl (E n+1 ) + T holc {E n+1 )) . The right hand side just defines the 
compact notation 0, and Q — > 1 in the classical limit. It is important to note that the Andreev states in this situation 
are degenerate. The state with electrons going to the right and holes going to the left has precisely the same energy 
as the state with electrons going to the left and holes going to the right. To be closer to our concrete model, let us 
introduce Tb, the travel time from one of the interfaces to the top, and Tt, the travel time through the top. The full 
traversal time for the case I is then 4Tf, + 2T t . 

In case II the electron coming from the left is reflected from the top, gets back to the left interface, undergoes 
Andreev reflection, gets to the right as a hole, gets back and finally undergoes yet another Andreev reflection. The 
energy spacing is again given by the inverse of the time of semiclassical motion. This time is now The Andreev 
states on the left side of the constriction do not overlap with the ones on the right side. This also leads to two-fold 
degeneracy of Andreev states for our symmetric setup. 

In the case III the electron coming from the left passes the constriction and undergoes Andreev reflection on the 
right side. The resulting hole goes to the left, is reflected from the top and undergoes Andreev reflection on the right 
side again. This results in an electron going to the left that undergoes Andreev reflection on the left side. There are 
thus four Andreev reflections involved. The state encompasses the electron and hole waves that propagate in both 
directions. Therefore it is not degenerate. The travelling time is now 8Tb + 2Tt- 

If T t <C Tb, the spacing of the states of the types I and II are evidently the same. For type III states, the travelling 
time is now approximately two times longer, which corresponds to a two times lesser energy spacing. However, there 
is no degeneracy in this case. This is why the number of states per energy interval remains approximately the same 
in all three cases. If T t 3> T b , there is a big spacing between type II states corresponding to energy scale En. The 
spacing of the type I and type II states is much smaller. 

Even a short glance at our numerical results proves that the picture based on the Eqs. (^|) and ( |To| ) is oversimplified. 
In the following section we present a more sophisticated picture, which accounts for the oscillating phenomena, found 
from the full solution of the Bogoliubov equations (0) . 



V. QUALITATIVE RESULTS 

Now we are in the position to present a concise summary of our qualitative results. Relevant parameters that 
determine the behaviour of the Andreev states are their energy E and the shift SE of the chemical potential with 
respect to the potential maximum. The superconducting gap A just determines the energy interval where Andreev 
states may be present, provided there are many states in this interval. 

The parameter regions exhibiting distinct behavior of Andreev states are presented in FIG. [| Two diagonal lines 
correspond to either electron or hole part of the state crossing the top of the potential relief. These lines separate the 
regions I, II and III discussed in the previous section. In these regions, the Andreev states are almost equally spaced, 
with the spacing being a slow function of energy. The energy scale of the spacing corresponds to min(Eii, En). The 
semiclassical theory relates the spacing with the classical traveling time. From this one would conjecture a sharp 
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crossover between the regions I, II, III, this presumably taking placeJn an energy interval of the order of several level 
spacings. This in fact corresponds to the situation described in Ref.El. 

The transitions between the regions I, II and III are much less trivial for a long shallow constriction described here. 
The reason for that is a strongreflection from the top of the constriction that takes place in a wide energy interval 
~ E c near the lines, see FIG. ||. The reflection takes place at two edges of the top. The interference between the 
reflected and transmitted waves disrupts a simple picture of semiclassical Andreev states. 

This is why we have three more regions, Ilia, la, and lb in the diagram. They are separated by dashed lines 
from their parent regions. Whereas the solid lines denote sharp transitions between the regions, the dashed lines 
correspond to smooth crossovers. For our concrete model, the crossovers are really very blurred, and the extra regions 
are difficult to separate from the parents. The reason for that is a slow decrease of the reflection coefficient with 
the energy, R w E c /E. The crossovers may become sharper in a different model of a long shallow constriction. The 
regions would remain the same. 

It is relatively easy to comprehend the situation in the region Ilia. There, the holes are completely reflected from 
the top and the transmitted and reflected waves do not interfere. They do interfere for electron components. This 
results in a regular modulation of spacings between Andreev states, which we call electron oscillations. We present in 
an Appendix a detailed theory of this effect. Here we give a final relation for "energy spacing times travelling time" 



G(E), see Eq. (A.13) 



o(e) = i + ^ cqs(2MVR-R (14) 

v ' T t + 2T b l + R-2cos(2(j) b )y/R V ; 

R being the energy-dependent reflection coefficient from the edge of the constriction middle and 4>b being a linear 
function of energy, <fib{E) = const. + ■nE/En. If Ej% <C En, each electron oscillation encompasses many Andreev 
states, of the order of En/E&. 

In the region la the situation is almost the opposite. Here, those are holes that are strongly reflected from the 
top. The electrons go at higher energy and are being transmitted much better. So one expects interference of holes 
only. We will call this phenomenon hole oscillations. Due to the slow decrease of the reflection coefficient mentioned 
above, the remaining interference of electrons is substantial. Therefore, the hole oscillations are not as regular as the 
electron ones. Also, the oscillation period is set by the inverse of the time for electrons to travel through the top, Tf" . 
The ratio of T t ft /Tf is not parametrically big, so that the oscillation encompasses only few Andreev states. In our 
numerical results, we see alternating spacing of Andreev states in the region la. 

In the region lb both electrons and holes exhibit reflections, and the interference picture thus becomes complicated 
and at best quasi-periodic. We will call this irregular oscillations. The visible quasi-period is En again. It arises as a 
result of interference of two reflections, one for the electron component and another for the hole one. These reflections 
are coupled by an Andreev process. One would not see such period in the interference picture of a normal electron. 

All these three types of oscillations are seen in our numerical results. Another effect of the reflection is that the 
degeneracy between Andreev states is lifted. For our symmetric setup, the Andreev states can be either even or odd, 
that is, symmetric or antisymmetric with respect to x — > —x. 
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VI. NUMERICAL RESULTS 



For numerics, we choose the transverse dimensions in such a way that only six modes can pass the constriction. For 
the sixth mode, the Fermi level matches the potential in the middle of the constriction and \8E\ -C H- We note that 
from now on A is chosen to be a real constant, since we assume no phase difference between the superconductors. If 
8E = only the electron can pass the constriction. As soon as SE < — A no transmission occurs for the waves forming 
Andreev states. Since the interesting phenomena occur for the sixth mode and for not too big values of 5E/A., we 
will show results for 8E/A taking values 10, 2, 1, 0.5, 0, -0.5, -1 and -2 only. 

Most calculations were done using a gap value of 4 x 10~ 5 /i, which is sufficiently small to achieve the semiclassical 
limit under investigation. 
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FIG. 6. O as defined by Eq. @ versus the relative Andreev bound state energy E/A, for a system length LkF = 7 x 10 . 
The first five modes support just one state each, while for the sixth mode the number of states is seven. Both electrons and 
holes traverse the system. 

In order to get a feeling for the system properties we first show the Andreev states spectrum for all modes and the 
two system lengths L = 7 x lO 4 /^ 1 and 2.1 x lO 5 ^ 1 in the FIGs. || and ^| respectively. For both cases 8E — 10A, so 
both components of an Andreev state can transit the system. The system lengths are such that the first five modes 
for L = 7 x lO 4 ^ 1 support one bound state each, while for the three times longer system two bound bound states 
are found. The energies of these states are plotted with the circles. 

The energy separation for the lower modes is approximately En, which is of the order of A choosen. The wave- 
functions of corresponding states easily pass the middle part of the constriction and their energies do not depend on 
8E, being very large for these states. Because they do not exhibit any interesting behavior we do not pay further 
attention to the first five modes. We just mention, that in all results to be displayed the Andreev approximation is 
made at the two SN interfaces. In the exact treatment the degeneracy of the two Andreev states with the two opposite 
propagation directions is lifted, and one ends up with two standing waves, a symmetric and an antisymmetric one. 
We calculated the exact energies for the five states depicted in FIG. ^| by filled circles, and found a splitting of the 
order of 10~ 10 /i, which is negligible indeed. 

More states are found for the sixth mode. Since they traverse the top, their energy separation is determined by 
Ei2 ^> En rather than En. We begin the discussion with FIG. (j| There are 20 Andreev states, as can be seen in 
the figure. Both the odd and even states are shown, the corresponding values being connected by broken and solid 
lines respectively. We see that the O values are close to 1, which proves the validity of the semiclassical treatment. 
So we concentrate on the deviations from 1. The energy distance from the top of the potential 8E « 47 E c . Although 
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the waves are still well transmitted, the reflection probability at the narrow part of the constriction becomes already 
appreciable (2%). 

This leads to two visible phenomena. First, the two states corresponding to two travelling waves propagating in 
opposite directions are mixed by reflection. By that the even and odd states have different energies. The stars in 
the figure indicate the result for an average energy, which corresponds to two degenerate states found after neglect of 
reflection at the constriction. The actual calculation in which only propagation to the right is allowed can be achieved 
by omitting the contribution of the wave propagating to the left in the normal metallic part at the left SN interface. 
The corresponding values are somewhat smaller than 1, but they lie at a stable height, in agreement with the 
estimate based on the forms (|9|), in which any reflection is not accounted for. The deviation of the stars from the 
lines is the measure of the effect. 
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FIG. 7. The same as in FIG. |(| but for a three times longer system length L/cf = 2.1 X 10 5 . The reflection has increased 
resulting in bigger amplitude of irregular oscillations of Andreev spacing. For the first five modes two states are found for each 
mode, while for the sixth mode the number of states is twenty. 

Second, the interference between the reflected waves leads to irregular dependence of on energy, as discussed in 
section The amplitude of these irregular oscillations is small corresponding to the small reflection coefficient. 

Both effects are better visible in FIG. ^ where the constriction is shorter and only seven states are seen. They are 
displayed by four triangles and three squares. However, this is not the main difference between the figures. The energy 
distance SE becomes w 23E C , which corresponds to a reflection coefficient of the order of 4%. The energy difference 
between the odd and even states becomes considerable, even such, that the highest state, a symmetric solution, does 
no longer have its antisymmetric counterpart. This would have got an energy E > A, while bound states are only 
found below the gap. The irregular oscillations arc increased in amplitude. 

Results for SE = 2A are shown in FIG. ||. In addition to the results for the system length L = 2.1 x 10 5 the even 
states are shown for the larger length of L = 7 x 10 5 . The number of states has increased to 119, and the oscillation 
frequency has increased as well. The two choices of L correspond to SE/E C w 10 and SE/E C w 21. The irregular 
oscillations are bigger in the first case, which corresponds to the bigger reflection. Due to larger number of states, 
one can see already a typical (quasi)period. The biggest interference scale En is 0.37A and 0.11A respectively. We 
see a good correspondence between this scale and the quasiperiod. 

The next figure, FIG. ^, presents the results for 5E = A, when Andreev states extend till the energy at which holes 
are reflected completely. We observe a new type of oscillations at energies > 0.6A. At these energies, the holes are 
reflected more effectively than electrons. So that the oscillations are determined mostly by their interference. The 
two choices of A correspond to SE/E C = 4.75 and SE/E C = 9.5. 
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FIG. 8. The same as in FIG. ^, but only for mode 6 states. Irregular oscillations in the region lb. The hole branch of 
the Andreev states lies always an energy higher than the gap energy above the barrier in the middle. We show results for two 
system lengths 2.1 x 10 5 and 7 x 10 5 , which give 38 and 119 bound states respectively. The reflection coefficient varies between 
0.05 and 0.1. The typical quaisiperiod is of the order of En. 



The most interesting figure is probably FIG. |l0|, which displays SE — 0.5A results. The parameter choice corre- 
sponds to SE ss 3.?>E C . All three types of oscillations can be simultaneously seen there. The oscillations are irregular 
till E pa 0.3A and they cross over to hole oscillations at higher energy. At E — 0.5A the holes are completely reflected, 
and the interference picture is determined solely by electrons. The transition is sharp indeed, the width being of the 
order E^ = 2 x 10~ 5 A. These electron oscillations are very regular. Their period is determined by En — 0.22A, each 
oscillation encompassing about 10 Andreev states. There is a striking contrast between the electron and hole oscilla- 
tions, the latter being much less regular. The point is that the hole oscillations are disrupted by residual interference 
of electrons and the electron oscillations are not, since holes do not penetrate the middle of the constriction and do 
not interfere. 

The electron oscillations throughout the full range of Andreev states are presented in FIG. [ll[ Here the width of 
the constriction is chosen in such a way, that U{x) in the middle coincides with the Fermi energy and SE = 0. All 
hole waves are reflected completely from the middle of the constriction, while all electron waves still can pass. The 
parameters correspond to A/E c ss 4.5 for the solid curve and A/E c 19 for the dashed one. This means that the 
reflection coefficient from the top varies from 1 to sa 0.2 or sa 0.05 respectively. This affects the oscillation amplitude. 
In FIG. |ll| one indeed sees that the position of the dips moves upwards with the energy, most clearly for the larger 
gap value. The spacing between the states and the oscillation period are determined by energy scales E^ and En 
respectively. This is why for four times larger gap four times more states and oscillations are found. 

We remind that up to now all calculations were done with the length of a flat part in the middle equal to 10% of 
the system length, L t = 0.1 L. We also show in FIG. [ll] the results for much shorter flat part length L t — 10~ 3 L. 
Six states are found, three odd states and three even states. For these parameters, the constriction middle is so short 
that Ei2 ^> En, and the spacing is of the order of En ps 0.37A. This is why the Andreev states for the short top 
match the oscillations of Andreev states for the long top. The matching is not perfect. The reason for that is the 
tunneling of the holes through the top, which becomes considerable for the short top (E^ = 0.5A). 

The following results are for parameter range SE < where holes are completely reflected from the top. FIG. |lj, 
with SE = — | A, is a counterpart of FIG.|o|, with SE = ^A. This means that the electrons are completely reflected 
at E < iA and transmitted at higher energies. By that the oscillations arise for the higher energies only, which is 
clearly seen in the figure. For E < ^A the two parts of the constriction are in fact not connected. The Andreev states 
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FIG. 9. For SE = A the bound states near A have a hole component touching the top of the potential. Strong reflection of 
these holes give rise to hole oscillations above « 0.7A. Two gap values correspond to 66 and 95 Andreev states. 

calculated are doubly degenerate, with one state resting in the left part of the constriction and the other state in the 
right part. The spacing between these states is of the order of En for both values of A. Again, this scale sets the 
oscillation period for the states above A/2, whereas their spacing is determined by much smaller energy scale E i2 . 
Two values of A correspond to A/E c « 4.5 and A/E c m 19. 

The results depicted in FIG. [f| correspond to SE — —A and SE = — 2A. So we have the space separation of 
Andreev states in the whole energy interval. No interference occurs, and the spacing is a smooth function of energy. 
The differences between SE = —A and SE = — 2A are small. For the smaller gap value the number of states is of the 
same order as the degenerate states for mode 6 in FIG. ^. The Andreev levels hardly depend on SE, since they never 
reach the middle of the constriction and do not experience the potential over there. 

VII. CONCLUSIONS AND PROSPECTS 

In conclusion, we present here a theory of Andreev states for a shallow long adiabatic constriction. The system 
differs from traditional models of S-N-S junctions by the absence of electron-hole symmetry and exhibits at least four 
disctict energy scales that determine the interference of the components of Andreev states. Such constrictions can be 
possibly realized in semiconductor heterostructures. 

We have started with a completely semiclassical approach to the states and reached in the first approximation a 
good agreemement with exact solutions of Bogoliubov equations. On the top of this, the exact Andreev levels may 
exhibit three distinct types of oscillating behavoir. These oscillations arise from interference of the transmitted and 
reflected waves in the middle of the constriction. We have extended the semiclassical theory to account for this 
interference and were able to reach quantitative agreement with numerical results. 

In this paper, we assumed no phase difference between the superconductors. It would be interesting to calculate 
the Andreev states in the presence of such phase difference and finally find the supercurrent in the structure. It is 
already clear that the supercurrent would exhibit an interesting dependence on controllable system parameters, in 
particular, on SE. This makes such calculation interesting in view of possible experiments. 
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FIG. 10. At 8E — 0.5A the hole component is completely reflected in the middle of the constriction at E > 0.5A. All three 
types of oscillations are seen simultaneously in this plot. 
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APPENDIX A: INTERFERENCE OF ANDREEV STATES 



The shape of the oscillations can be understood in quasiclassical terms if the interference of the waves is taken into 



account. The configuration of the system and the parameters used are shown in FIG. 14 



To evaluate the quantization condition, we consider the amplitude of the wave that starts propagating from point 
X, goes all the way through the constriction and gets back to the same point. The starting point can be chosen 
arbitrarily. Each part of the trajectory contributes a factor of the form e*^ to the amplitude. The phases (j>i, 02 and 
04 are shown in the figure, 03 is the propagation phase from B to C, and 05 is the phase from D to X. If there is no 
reflection at A, B, C and D 



Amplitude = e^e^e^e^e** 6 



(Al) 



in which the second equality is the condition for an Andreev bound state. This quantization condition can be written 

as 



+ d>4 + d>5 = 2wn 



Using the symmetries 



h + 05 = 03 = 206 and <j) 2 = 



>4 = <pt 



(A2) 



(A3) 



one can write down the equality 
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FIG. 11. The electron oscillations through the whole energy range. The hole component is completely reflected for all 
energies. Each oscillation encompasses about ten states. Circles and stars represent the states for a very short middle of the 
constriction and for the smaller gap value of A = 4 x 10 _5 /i. Note the correspondence of their positions and the oscillation 
dips. 
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2T t + AT b , 



(A4) 



where T t is the time for the classical motion between the points A and B, and Tf, is the time it takes to get from the 
interface to the constriction top. Here we made use of the quasiclassical relation between the phase and the time of 
the classical motion = T. Eq. ( |A4| ) can be considered as a slight generalization of Eq. (|l0|). 

Now we extend the analysis by assuming the possibility of relections at the points A, B, C and D. For each point 
four coefficients are to be defined, namely tin and t^L for transmission from left to right and from right to left 
respectively, and r^n and trl, similarly, for reflection. By summing up all possible ways of propagation, from the 
point X to the right and back to that point coming from the left, and accounting for all reflections and transmissions 
one arrives at 



A.RtftT.A-AD 



Amplitude = e'* +fc >[rS + 



l-r^A AD 



= 1, 



(A5) 



in which the second equality is the quantization condition. The quantity A^ is the amplitude to get from A to D 
and is given by 



,(2) .(2) i0 3 

A A D = e«* a+ <*')[r™ + tL ^ e 



1 - r^e^3 ' 



Again using the phases 4>b and 4>t, solving Eq. (|A5|) for A ad, and using Eq. 



a 2i<t> t 



r (l)„i206 



LR 



O i20(, 7 ,( 1 ) 



' LR 



(A6) 

), one obtains the symmetric form 

(A7) 



RL 



in which and are defined by 



13 



0(E) 



1.3 



1.2 



1.1 



1.0 



0.9 



0.8 



0.7 



0.6 



0.0 



5E/A=-0.5 
Lk F = 3.5x1 5 



X X X X X X 



x 



0.2 



*A/u.= 4x10 
- A/n = 4x1 0" 5 
x A/u, = 1.6x10" 
A/ji = 1.6x10" 




0.4 0.6 

E/|A| 



0.8 



1.0 



FIG. 12. In this case, E — 0.5A separates regions II and III b. The Andreev states in the region II are localized near the 
superconductors and do not go through the middle of the constriction. Note the correspondence between the spacing of the 
states at low energy and oscillation period at higher energy. 



_(l,2) .(1,2) (1,2) (1,2) (1,2) 
* — LR RL ' LR 1 RL ' 



(A8) 



We now study the general expression (A7) using simplifying assumptions. If we disregard scattering phases and their 
energy dependence we can write 



r LR. ~ r LR ~ r RL ~ r LR. 



R, 



(A9) 



R being the reflection coefficient. By that we find, that z^ 1,2 ) = —1 and the quantization condition (A7) becomes 

r >, ..... j \ /?r 



1 - e a <t>»y/R 



which is equivalent to 

2(4> t + 20 6 ) + 4 atan 
For the energy derivative this gives 

dn 



sin(20 b )Vi? 



2nn. 



1 - cos(20 b )Vi? / 

cos(2(j) b )y/R - R 



2nn— = 2T t + AT h + 4T b - _ . 

dE 1 + R-2cos(2<f> b )VR 

We rewrite this using our definition of and obtain the desired expression, see Eq. (|l^), 

e(£) = 1 , 2T b cos(20 fc )Vi?-i?^ 



T t + 2T b l + R -2cos(2<f> b )VR 



(A10) 



(All) 



(A12) 



(A13) 



If there are no reflections, 9 is equal to 1. If reflections are accounted for, Eq. ( Al3| ) implies, that for R <C 1 small 
modulations below the value of 1 are found in plotting the quantity 5E(T t + 2Tb)/n as a function of E or E/A. If 
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FIG. 13. Localized Andreev states in the region II. Their spacing is a smooth function of energy. Their energies hardly 
depend on the position of the chemical potential with respect to the maximum of the potential. 



R ~ 1 periodic dips below the value of 1 are coming out, completely in agreement of the structure of FIG. |llj. For 
a constant R the dips lie at equal height. In reality R is energy dependent at E ~ E c and decreases with increasing 
energy. Recalling the definition of 4>b,4>t, we see that the spacing is determined by the energy scale E.^ (provided 
E ~ E c ) and the oscillation period is determined by the bigger energy scale En. 

Similar consideration can be performed for hole oscillations in the region la. This results in oscillations with a 
period determined by difference of phases accumulated respectively by electrons and holes during their travel though 
the middle of the constriction. 

The irregular oscillations in the region lb can be also quantified along these lines. However, here the reflection 
coefficients for the electron and hole wave are different, and different phase differences are accumulated by electrons 
and holes in the constriction middle and beyond. The typical quasiperiod is determined by the phase <frb, that varies 
with energy in a slowest way. 
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FIG. 14. Schematical picture of our constriction, including the paths of the electron and hole branches of the Andreev bound 
states. The Fermi energy coincides with the flat part of the potential at the constriction. In the quasiclassical derivation 
reflections of the electron wave are accounted for which are taken to occur at the the entrance and exit of the constriction. 
Reflection 1 occurs at the point A, which coincides spatially with point D, and reflection 2 occurs at point B (or C). 
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